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I. INTRODUCTION 



Generalized parton distributions (GPDs) [1-3] are now the object of intensive theoretical studies, especially within 
the context of applications to deeply virtual [4,5] (for a detailed recent review see Ref. [6]) and large momentum transfer 
processes [7,8]. The main advantage of GPDs is their universality allowing to connect different hard processes, both 
exclusive and inclusive. The price for this is the complexity of GPDs: they are functions of 3 variables, e.g., skewed 
parton distributions (SPDs) J^i^{X,t) or H{x,^;t) depend on the fraction X (or x) of the momentum carried by the 
active quark, the skewedness parameter ( (or ^) and the invariant momentum transfer t. For this reason, the most 
promising approach to disentangling GPDs from experimental data is to construct realistic models for GPDs and 
fix their parameters by fitting the data. The crucial point for the model building is that, in specific limits, GPDs 
reduce to more familiar functions describing the hadronic structure, such as usual parton densities, form factors and 
distribution amplitudes. The "reduction" relations between GPDs and these functions have been used as a basis for 
building phenomenological models of GPDs [9] . Another fruitful idea used in the model- building is to construct GPDs 
from the light-cone (LC) wave functions [7,8]. The most popular Ansatz [10] assumes a Gaussian dependence of the LC 
wave functions ip{x,k±) on the transverse momentum k±. A pragmatic reason behind this choice is the simplicity of 
Gaussian integrals allowing to obtain many results in analytic form. However, there are no a priori grounds to exclude 
wave functions with other types of transverse momentum dependence. In particular, the two-body {i.e., qq) component 
of the pion wave function was calculated recently in a model [11] based on the one-gluon exchange approximation in 
the light-front framework. The wave function was found numerically, and it was observed that the fit is better if one 
uses a power-law form rather than a Gaussian [11,12]. Furthermore, the power-law wave functions were used some 
time ago in models for the nucleon form factors [13]. In the present paper, we show that a simple power-law Ansatz for 
the pion LC wave function allows one to obtain explicit analytic expressions for the form factor and generalized parton 
distributions. To make our presentation self-contained, in Section II we remind basic information about generalized 
parton distributions which is used in the following Sections. In Section III, as a starting example, we consider the 
model with the Gaussian dependence of the LC wave functions on the transverse momentum. In Section IV, we 
specify the explicit "toy" model expression for the effective pion wave function, which is then used in Section V to 
derive a parametric representation for the pion form factor. We show that the two parameters of this simple model, 
the constituent quark mass and the wave function width, can be easily adjusted to provide a curve close to existing 
experimental data. In Section VI, we analyze the asymptotic large-Q^ behavior of the pion form factor. We consider 
both the massive m ^ and massless m — cases. We show that, in the latter case, the pion form factor in our model 
with a power-law wave function il}{x,ki_) ~ l/[-\/x(l — x){l + k\/{}?x{l — x))"] has the same asymptotic behavior 
Ftt{Q'^) ~ l/<3^ for any power n. Though this behavior is generated by the soft (Feynman) mechanism, it formally 
coincides with the quark counting law dictated by the hard one-gluon exchange mechanism. We show that the 1 / 
behavior of the soft contribution is related to the fact that the parton distribution f{x) in the massless case does 
not vanish for x = 1. In massive case, f{x) — > as a; — > 1 and the ultimate asymptotic behavior is ln((5^/A^)/(5^. 
However, for a wide range of accessible Q^, the curve mimics the behavior. In Section VII, we note that our 

parametric representation for the form factor has the form of the reduction relation connecting the pion form factor and 
the double distribution (DD) F{x, y; t). The DD obtained in this way has correct spectral and symmetry properties. 
Moreover, it has the factorized structure proposed in Ref. [9] : it looks like a distribution amplitude with respect to the 
y variable and like a parton density with respect to the x variable. It also provides a nontrivial example of the interplay 
between x, y and t dependence of DDs. With an explicit model for DDs at hand, one can calculate the relevant skewed 
distributions: the nonforward parton distribution T(^{X; t) or Ji's off-forward parton distribution H[x, ^, t), see Section 
VIII. In the simple toy model that we use the pion is treated as an effectively two-body system, which is not very 
realistic: one may expect that the parton densities at small x are affected by many-body components. Indeed, the 
valence parton density obtained in our model differs rather strongly from the phenomenologically established form. In 
Section IX, we propose to fix this deficiency by adopting a model with a more realistic x-profile at t = 0, but preserving 
the analytic structure of the interplay between x, y and t dependence generated by the power-law Ansatz. We show 
that by slightly changing the quark mass and the wave function width parameter it is still possible to get a good 
description of the pion form factor data. We present SPDs J^(^{X;t) obtained from the "realistic" DD. In particular. 



2 



we show that in the "soft pion hmit", C = l,t = 0, the isovector part of the "reahstic" SPD has the shape close to 
the asymptotic form of the pion distribution amphtude. In Appendix A, to demonstrate that the variables x, y of the 
parametric representation for the form factor indeed have the meaning of the variables of double distributions, we give 
a covariant derivation of the toy model DD in a scalar model. In Appendix B, we discuss the structure of model SPDs 
in the impact parameter representation. In particular, we show how one can use superpositions of power-law DDs to 
build models for SPDs satisfying positivity bounds. In Appendix C, using again the toy scalar model, we briefly show 
how one can use our approach to build the models for two-pion distribution amplitudes that appear in 7*7 tttt 
reaction, which can be treated as the crossed-channel process to deeply virtual Compton scattering. Our conclusions 
are formulated in Section X. 

Summarizing, in this paper wc construct power-law models of the C-odd double distributions F(x, y, t) for the pion 
and the relevant skewed parton distributions J^^{X;t). By construction, the model GPDs satisfy such important con- 
straints as reduction relations to usual parton densities and form factors, they have correct spectral and polynomiality 
properties, thus providing a model that can be used in phenomenological applications. For the simplified scalar case, 
we also build the models that automatically satisfy the positivity constraints. 

II. BASICS OF GENERALIZED PARTON DISTRIBUTIONS 

Generalized parton distributions parametrize nonforward matrix elements of composite operators. To define the 
leading-twist GPDs for the pion, we start with 

i-{P - r/2 I D,,...D,Ji^a\P + r/2) = 2 ^ 4fe (t) {P,P,. ■ ■ ■ V.-.+i • • • V J 

+l-D(-)it){r,r,,...r,J , (2.1) 

where D— {D — D)/^, {. . .} denotes the symmetric-traceless part of a tensor, a numerates quark flavors and the 
quark fields are taken at the origin. Compared to the more familiar case of forward matrix elements defining the usual 
parton densities, we have two 4-vectors, P and r, both of which can be used to build the tensor structure of the right 
hand side of Eq. (2.1). The index k specifics how many times the vector r appears in a particular term of the sum. 
Incorporating hermiticity properties of the local operators and time-reversal invariance, one can show [14] that k is 
even. Now one can define double distributions f{p,a;t) as functions generating Ank{t) through its f3"'~''a'' moments 

{l±(-l)"}4«Hi)= / d(3 p->'a'f^{P,a;t)da. (2.2) 

J-l J-l+\f3\ 

The spectral property -|- |q;| < 1 can be proved for any relevant diagram of perturbation theory [1,5]. 

As usual, the Mellin moments define two functions: f~{/3,a;t) corresponds to even n while f^{0,a;t) to odd n. 
They both are even functions of a. With respect to P, f^{P,a:,f?j is even while f^{P,a-,t) is odd. For /3 > 0, one 
can write f~{P,a;t) as the difference fa{P,a;t) — fa{P,a;t) of quark and antiquark distributions {i.e., f~{P,a;t) 
corresponds to a valence quark distribution: /~ = ) and f^{P,a;t) as their sum fa{P,a;t) + fa{P,cx;t). The 
Polyakov- Weiss D-term [15] is defined as the function Da{a;t) whose a" moments give the Di^\t) coefficients. The 
latter are nonzero only for odd n, hence Da{a;t) is an odd function of a. 

We stress that this definition of double distributions is absolutely Lorentz invariant: it does not require reference 
to any particular frame. Moreover, the mutual orientation and relative size of two momenta P and r are arbitrary. 
If, in some particular frame, the space part of the momentum P is oriented in the (longitudinal) a;3-direction, the 4- 
momentum r may also have a nonzero longitudinal component, but it may be purely transverse as well, having nonzero 
components in the transverse a;i,X2-plane only. The double distributions f{P,a;t) parametrizing the nonforward 
matrix element are Lorentz invariant objects and they are the same in all cases. 

Usually, to extract the symmetric-traceless part of a tensor 0^^j...^„, it is multiplied by z^zf^^ . . . z^'^, where z^ is 
a lightlike vector z^ = 0. This trick corresponds to a projection of Eq. (2.1): 
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{P - r/2 I i,,z{iz B)"Va I P + r/2) = 2{Pz) ^ ^"^(i) (P^)"-'=(r^)'= + -i)(«)(f)(r^)"+i , (2.3) 

fc=0 

(where z = jfj.z'^). The direction of z is arbitrary, but, to access all the coefficients Ank{t), one should have both 
(Pz) 7^ and (rz) ^ 0. In particular, if z has only the minus light-cone component, both P+ and r+ should be nonzero 
to make all the coefHcients visible. Such a situation is characteristic for deeply virtual Compton scattering (DVCS) 
where the momentum transfer r must have a nonzero longitudinal component. To study DVCS, it is convenient to 
treat the ratio ^ = r"'"/2P+ as an independent variable and define off-forward parton distributions H{x,^;t) [2]. To 
this end, one introduces the functions 

n I 

= i: k\{n:ky Ai^^it)^'+Di:\t)e^' (2.4) 

and declares M.l^\^; t) to be the moments of Ha{x, ^; t): 

{l±{-ir}Mi^\tt) = j\^H^{x,^-i)dx . (2.5) 
These definitions provide a formal relation between H{x, ^; t) and /(/?, a; t) 

/•I r^-m 

Ht{x,^;t)= dp f^{P,a;t)5{x -P-^a) da +{l±l)sign{0 D,{i:/i;t). (2.6) 

J-l J-l+\l3\ 

Combining Eqs.(2.2) and (2.3) gives the definition of DDs through the parametrization of nonforward matrix ele- 
ments of nonlocal light cone operators 

{P-r/2\M-^/2)zM^/2)\P + r/2)U.^o = {Pz) / dp e-'^^''^^-"'^^^^'^ (f+ {P,a;t) + f- {P,a;t)) da 

J-i J-i+\0\ ^ ^ 

+ {rz)j e-'"(''^)/2Z)„(a;i)da . (2.7) 

Using the symmetry of ftiP-, o.\ t) with respect to P and a, one can rewrite this representation in terms of quark 
and antiquark DDs taken for positive P only 

.1 .1-/3 

(P-r/2|^„(-0/2)5^„(^/2)|P + r/2)U2=o = 2(Pz) / dp 

Jo J-i+p 

X (/„(/3,a;t)e-*'^(-P^)-'«('-^)/2 - /<i(/3, a; t)e*'5(-P^)+'«(")/2) da 

+ {rz)j e-*«('-^)/2£)„(Q!;i)rfa . (2.8) 

In the forward limit, r = 0, the left hand side of Eq. (2.8) coincides with the matrix element defining the usual 
parton densities fa,a{x)- This gives the reduction relations 

/I— X 
fa,a{x, a;t = 0)da = fa,a{x) and i?a,a(a;, ^ = 0; i = 0) = fa,a{x) . (2.9) 
-l+x 

On the other hand, keeping r ^ but taking n = in Eq. (2.1) one deals with the matrix element of the vector 
current which defines the a-component Fa{t) of the relevant form factor. The reduction relations connecting GPDs 
with form factors result from 

P„(i)=4o^i) (2.10) 
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and are given by expressions 



^ '^(^UP,a;t)-Mf},a;t))da = F,{t) , ^ (^H,{x,^;t) - H^{x,^;t)) dx = Fait) (2.11) 

containing only the valence quark combinations — fa ^ fa and H™^ = Ha — Ha- 

The representation (2.8) has the structure of a plane wave decomposition, which provides the parton interpretation 
of DDs: the quarks carry the momentum /3P + (1 + a)r/2 originating both from the average momentum P and 
the momentum transfer r. Another possibility (which is more convenient in applications involving light-cone wave 
functions) is to write the momenta of quarks as xpi + yr, i.e., in terms of r and the original hadron momentum 
Pi = P + r/2. The new variables x,y are expressed through P, a hy x = l3,y = (1 + a — /3)/2. The resulting DDs 
Fa,a{x,y;t) "live" on the triangle < x,y,x + y < 1. Since f{f3,a;t) are even functions of a, the DDs F{x,y;t) 
are symmetric with respect toy^ 1 — x — y transformation ( "Munich" symmetry [16]). For light-cone dominated 
processes, like DVCS, only the plus component xp^ +yr~^ is essential. Defining the skewedness parameter ( = r+/p^, 
we introduce nonforward parton distributions [5] 

Tl^~\X,t)= j^dx F''^%x,y;t)5{X-x-Cy)dy . (2.12) 

These distributions are related to usual parton densities by 

FaAx,y;t = Q)dy = fara{x) , Tl'^o{^;t = Q) = f^A^) (2.13) 



/ 

Jo 
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and to form factors by 



[ dx [ F™' (x, y;t)dy = Fa [t) , / [X, t) dX = F„ (i) . (2.14) 

Jo Jo Jo 

Note that the double distributions F{x, y; t) are integrated over y in both of the above reduction relations. Thus, it 
makes sense to introduce intermediate functions 



^a,a{x,t)=l Fa,a{x,y;t)dy = J^c^^o{x;t) . (2.15) 
Jo 

They satisfy simpler reduction relations 

:Fa,a{x, i = 0) = fa-a{x) and / ^^^2;, t) dx = Fait) . (2.16) 

Jo 

Thus, the functions J^ix,t) are hybrids of form factors Fit) and usual parton densities fix), that is why we call them 
nonforward parton densities (NDs) [7] . 



III. GAUSSIAN WAVE FUNCTION AND NONFORWARD PARTON DENSITIES 

The concept of NDs can be easily illustrated within the framework of the light-cone formalism. Consider a two-body 
bound state whose lowest Fock component is described by a light-cone wave function ipix,k±). In a frame where the 
momentum transfer r is purely transverse r = r±, one can write the two-body contribution into the form factor as [17] 

F^'^^\t)= [ dx [ ^*ix,k_L + il-x)r±)ijix,k±)d'^k±= [ T^'^^\x,t)dx , (3.1) 
Jo J Jo 

where F^'^''\x,t) is the 2-body contribution into the nonforward parton density 
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jr(2b)(a;,t) = / tjj*{x,k± + {l-x)r±)ip{x,k±)(fk± 



(3.2) 



Adding contributions from higher Fock components, one obtains the total ND J^{x, t) whose integral over x gives the 

form factor F{t) of the bomid state. As discussed in the previous section, at zero momentum transfer T{x,f?) reduces 
to the usual valence parton density f{x) = !F{x;t = 0). Furthermore, there is the usual form factor normalization 
condition F{t = 0) = 1. Finally, for the valence quark distributions, the integral of f{x) over x is 1. These conditions 
are satisfied in the simplest way by the factorized Ansatz T{xJ,) = f{x)F{t), in which there is no interplay between 
X and t dependence of T{x,t). One may expect that in reality the situation is more complicated. Consider a wave 
function with a Gaussian dependence on the transverse momentum k± (cf. [10]) 

^{x, k^) = (^(a;)e-'=i/"^(i-^)^' (3.3) 

(note, that k'^/'lx{l — x) is essentially P written in the light-cone variables x, k±). Taking the Gaussian integral over 

x,k± we get 

JP^^b) (a;, t) = (x)e(i-^)*/4^^' , (3.4) 

where 

/(2'')(x) = 7rx(l - x) h? tp'^ix) = :rW(a;, t = 0) (3.5) 

is the two-body part of the relevant parton density f{x). To get the total result for either usual f{x) or nonforward 
parton densities T{x,t), one should add the contributions due to higher Fock components. These contributions are 
not small, e.g., with the Gaussian Ansatz t]ic valence du contribution into the normalization of the tt"*" form factor for 
f = is about 25% [10]. The problem is that we do not have a formalism providing explicit expressions for an infinite 
tower of light-cone wave functions. However, the parton densities f{x) are known from experiment. In this situation, 
one can treat Eq.(3.4) as a guide for fixing interplay between x and t dependence of NDs and model them by 

J-«(a;,t) = /„(x)e^*/^-^\ (3.6) 

The functions fa (x) here are the usual valence a-quark parton densities. One can take them from existing parametriza- 
tions of parton densities like GRV, MRS, CTEQ, etc. This model (originally proposed in Ref. [18]) was succesfuUy 
applied in Ref. [7] to describe the proton form factor Fi{t) in a wide region 1 < — t < lOGeV^ of momentum transfer 
by fitting the only parameter A^ characterizing the effective proton size. 



IV. POWER-LAW WAVE FUNCTIONS 



GPDs give the most general parametrization of nonforward matrix elements. Furthermore, both of them, DDs 
F{x,y;t) and SPDs J^i^{X;t) are functions of 3 variables: in addition to the invariant momentum transfer t they 
depend on two "longitudinal" variables x, y or A, C^. However, the Gaussian model of the previous section gives a 
representation for the form factor in terms of a one-dimensional x-integral of the function J^{x, t) depending on only 2 
variables, x and t. One may suspect that the Gaussian Ansatz is a degenerate case failing to reveal the richer structure 
present in more general situations. In what follows, our goal is to study a model based on power-law wave functions. 
As we will see, though this model is more complicated, we are still able to get most of the results in analytic form, 
which allows us to use it for building nontrivial Ansdtze for generalized parton distributions. 

The qq wave function of the pion found numerically in [11] was parametrized analytically by a power law fit 




with K ~ 2 rather than by a Gaussian ip{k) ~ exp(— fc^/A^). Here k"^ = k"^ + k\ is the square of the relative 3- 
momentum and A is the parameter characterizing the width of the A;^-distribution, i.e., the size of the system. In 
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the case of equal quark masses, there is a simple relation [19] between the usual variables {kz,k±) and the infinite 
momentum frame (IMF) variables {x,k±) 



1 

^=2 



(4.2) 



where m is the effective quark mass. The relation between ip{k) and the IMF wave function il>{x, k±) is given by 

^{x,k^) = ipik) ' -. (4.3) 

Va;(l - X) 

For light quarks, one may expect that the size parameter A is close to the effective quark mass m. Then the factor 
(l + Zc^/m^)^/^ can be essentially absorbed into a redefinition of the power k, whose precise value, in fact, is not critical 
for our purposes. Thus, in what follows, we will consider a simplified power-law IMF wave function 

ipix, k_i_) = ^ , (4.4) 

' y'x{l-x)[a + bkl]^ ^ ' 

where a + bk^ is the IMF version of (1 + k"^ / h?) with 

^_ ^j^s^i^- - \Y ■^^ s-^ (4 5) 

x{l — x) ' Am?x{l — x) ' A ' 

and N is the normalization constant. For the 2-body Fock component, it can be fixed from the requirement that the 
integral of tl}{x, k±) over x and k± should give the pion decay constant 



/ V'^^'Ha;, k^)dxdPki_ = U ■ (4.6) 

As noted in the previous section, the knowledge of the 2-body wave function is not sufficient to calculate the pion 
form factor. To get it, we should add the contribution from all higher Fock components. Just like in the case of the 
Gaussian wave function, the two-body component is responsible only for some portion of "1" in the normalization 
condition -F(O) = 1 [10,12]. Again, the structure of higher Fock components can be only guessed. To avoid making 
too many guesses, we will analyze the simplest "one-guess" model in which a single two-body-like function tjj{x,k±) 
(4.4) imitates the contribution of all Fock components into the pion form factor. Thus, we take 

F{Q^) = j i^{x,k^ + {l-x)q^)tlj{x, k±)dxd^k± , (4.7) 

and normalize the effective wave function tp{x,k±) by 

^{x,kjL) \^ dxd'^k± = l . (4.8) 

This condition (4.8) gives an explicit expression for the normalization constant A'' of the effective wave function: 

where 

./X dx (1-^2)3 

Before proceeding further, we would like to make it clear that substituting the total contribution of higher Fock 
components by a 2-body type term is just a toy model, and we do not expect it to adequately describe all the aspects 
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of the pion structure. In particular, the total parton density in the toy model has the same {x ^ 1 — x symmetric) 
shape as its 2-body part, and it vanishes at x = 0. One would expect, however, that the contributions of higher Fock 
components are shifted to smaller and smaller x values producing eventually the experimentally observed ~ i/\/x 
behavior. We do not know how much each term of the infinite tower of Fock components contributes to the parton 
density, but we know (from experiment) what is the total result. Thus, our ultimate strategy, just like in the case 
of the Gaussian wave function, is to calculate GPDs in the toy model, identify the factor corresponding to the usual 
parton density and substitute it by the experimental one. On the other hand, one may expect that form factor, being 
an integral of the relevant GPD, should not be too sensitive to the details of its a;-dependence, at least in some range 
of momentum transfer t. Thus, we study first the form factor in our toy model. We show that, despite its crudeness, 
the toy model can easily fit the form factor data by adjusting the two parameters of the model. Then wc incorporate 
the main advantage of the toy model, the possibility to do calculations analytically, and obtain the representation for 
the form factor in terms of DDs. Finally, we "correct" the latter in such a way that, after integration, they produce 
experimental parton densities. We also show that this model gives DDs with a nontrivial "profile" dependence on the 
y variable. 



V. FORM FACTOR IN TOY MODEL 

The k± integral in the expression for the form factor 



F(O^) =N^ f ^-^ (5 1) 

' J x{l-x)[a + b{k^ + {l-x)q^ma + bkl]^ ^ ' 

can be done using either the Feynman parameters or the Schwinger a-representation method briefly described below. 
To this end, we use 
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A- TiK) Jo 



a"-^ e-""^ da , (5.2) 



where k = 2 in our case. After calculating the Gaussian integral over k±, we arrive at the representation for the form 
factor in terms of two parameters ai and a2 

^ ' Jo Jo ("1 + "2) Jo bx{^ -x) ' ^ ' 
where Q"^ = Q±- Changing the variables 

ai+ a2 = A, ai = 7A, 0:2 = (1 — 7)A, daida2 = XdXdj , (5.4) 
we obtain the parametric representation 

F(Q2) = [' rf^^(i _ ^) r a2 dAe-I^^+^^^i-^^'^^i-^^^'l . (5.5) 

Jo bx(l - X) Jo Jo 

Integration over A is easily performed to give 

Incorporating the normalization condition, we get the final result 

FAQ') = f'dx r d7 f ^n213 • (5-7) 

Ms) Jo Jo [a + 6(l-x)27(l-7)Q2]3 

By construction, the form factor has the correct value at = 0. However, its slope at this point depends on the 
values of the model parameters m and s. To obtain the analytic expression for the slope we note that, in the small 
(52 limit, one can expand the denominator of the 7 integration 
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FAQ') Iq^^o =4tJ''^ f ^77(1 -7(1 - l)Q\l " xf 



Ms) Jo a3 7o 
6&' 



+ ^(7(i-7))V(i-:rr + 



(5.8) 



Using the normalization condition (4.8) for the wave function, taking the integrals over 7 and introducing the variable 
z = 2x — 1 (like in Eq. (4.10)) we obtain 



where 



and 



FAQ': 



I ^ 3 .Bjs) 9 

20 m2 A(s) 280m4 y4(s) 



1 (l + z2)(l_;22)3 



Cis) 



S2)z2)4 
1 (l+z4 + 6z2)(l-22)3 



(1-(1-S2)z2)5 



dz 



(5.9) 



(5.10) 



(5.11) 




).2 0.4 0.6 0.8 1 
FIG. 1. Combinations s'^ B{s)/A{s) (solid) and s'^C{s)/A{s) (dashed) as functions of the parameter s. 

The integrals A{s), B{s) and C(s) can be calculated in elementary functions, though the results are rather lengthy. 
Fig. 1 shows the plot of the combinations s'^B{s)/A{s) and s^C{s)/A{s) as functions of the parameter s. The 
combination s'^B{s)/A{s) is monotonically increasing from zero to infinity. Hence, after choosing the effective mass 
m we can always find such a parameter s that the slope dFT^{Q^) / dQ^ of the pion form factor at Q2 = q has the 
experimental value dF^^^ [Q"^) / dQ"^ ~ 1/™p [20]. For masses m — 0.2 , 0.3 and 0.4 GeV, the parameters s fitting the 
pion charge radius are 0.56, 0.95 and 1.33, respectively. 




0.2 0.4 



0.6 
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FIG. 2. Left: Form factor -Pn-(Q^) for three different parameterizations of the wave function: with m = 0.2 GeV (solid), 
0.3 GeV (dashed) and 0.4 GeV (dotted). Right: Q'^F^iQ'^), with the same wave functions as in the left. 

In Fig. 2, left, we have plotted the form factor as a function of in the low region < 1 GeV^ for these three 
different parametrizations of the wave function. Since they have the same slope at Q"^ = 0, the curves are rather close 
to each other. However, the difference between the curves becomes more pronounced as increases. In Fig. 2, right, 
the form factor calculated with these three effective wave functions is shown in the region up to 5 GeV^ relevant 
to future experiments at Jefferson Lab. Among these three choices, the most close to existing experimental data [21] 
is the curve corresponding to m = 0.3 GeV and s = 0.95. 



VI. ASYMPTOTIC BEHAVIOR OF THE PION FORM FACTOR 

According to Fig. 2, right, in the accessible energy range < 5GeV^, the model curves show the behavior close 
to the 1/Q^ scaling. Since the mass scales involved are rather small: m^, A^ 0.1 GeV^, one may think that the form 

factor is already in the asymptotic region. 

To analyze the asymptotic behavior of the form factor one can follow the approach described in Ref. [22] . The basic 
idea is that in the Drell-Yan formula (4.7) we deal with an overlap of two functions ■ip{x, k±) and ip{x, A;_l — (1 — x)q±) 
whose k±_ arguments arc separated by a gap (1 — x)Q in magnitude. Furthermore, ip{x,k±) rapidly decreases with 
k± . Hence, when is large, the integral over k± in the form factor expression is dominated by two regions of phase 
space [22]: 

1) I kj_ [<§; (1 — x)Q, where ^jj{x, k^) is large; 

2) I A;_L + (1 — x)q± ]<C (1 — x)Q, where tjj{x, k± + {1 — x)q±) is large. 

In the first case, k± can be neglected compared to (1 — x)q± in the wave function. The contribution from this region 
is then approximated by 

Fl^\Q^)r^J^ i;{x,{l-x)Q)dx J 0{\k±\<{l-x)Q)^{x,k±)d^k^. (6.1) 

Since the wave function falls off rapidly at large transverse momenta, the major contribution to the integral comes 
from the region where | A:^ | is much smaller than (1 — x)Q, and one may hope that the k± integral of ipi^i k±) can be 
approximated by the pion distribution amplitude 'p{x). The next statement usually made is that the large-Q behavior 
of the function tjj{x, (1 — x)Q) is determined by the large-fcx behavior of ip{x,k±_) and, hence, the large-Q behavior 
of the form factor just repeats the largc-fcj^ behavior of the wave function; i.e., if tp{x,k±) ^ 1/fc", then F ^ l/Q"". 
Note, that the last statement is only true if, after these substitutions, the integral over x converges. However, it is 
easy to derive that, after the k± integration in Eq.(6.1), the remaining integrand for the a;-integration is proportional 
to 



c{l-x) + s^ {x-lf + {l-xr^' 

Neglecting the x{l — x) and s'^{x — 1/2)^ terms compared to (1 — x)'^Q'^ /AK^, one would get the integral dx/{l — x) 
logarithmically diverging in the a; — > 1 region. Of course, this approximation is only true when (1 — x)'^Q'^ /Ah? » 1 
or a; <C 1 — 2K/Q. This cut-off converts the logarithmic divergence into ln((5^/A^). Hence, the asymptotic behavior is 
~ ln((5^/A^)/Q^. This result can also be obtained from our representation (5.7) for the form factor, which we write 
now as 

1 /■ ^ /■ "^"-^"''"-''>° . (6.3) 

Ji(s)Jo Jo [^(1 - I) + - ipr- + t(1 - 7) "~4Ay ] 

Again, there are two regions: 7 ^ 4A^/[(1 — x)'^Q'^] and 1 — 7 <C 4A^/[(1 — x)^(5^] producing the leading large-Q^ 
contribution ~ A^/[(l — x)'^Q'^]'^. Combined with other x-dependent factors, this gives the dx/{l — x) divergence 
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or, after a more accurate calculation, the \u.{Q'^ / A?) / Q'^ asymptotic behavior. This behavior is not yet visible in the 
curves shown in Fig. 2, right. The curves suggest, in fact, the behavior. The slow approach to asymptopia can 
be traced to rather small numerical factor 7(1 — 7)/4 ^ 1/16 accompanying the term. As a result, the effective 
scale governing the Q^-behavior is something like 16A^ ~ 1.5 GeV^ rather than simply A^. Thus, the quark mass 
squared rm? ^ O.lGeV^ is small compared to the effective scale, and it is worth investigating what happens when 
quarks are massless, i.e., when s = 0. Then 



67(1 -7)x^ 



-7(1-7)^ 



3 ' 



(6.4) 



and situation drastically changes: the large-Q^ behavior is dominated by integration over 1 — x <C 4A^/[7(1 — 7)Q^] 
region. The remaining integral over 7 has no singularities, and we get ^^(Q^) ~ A^/Q^ for the asymptotic behavior. 
Clearly, the asymptotic behavior for massless quarks is governed by the soft (or Fcynman) mechanism. One can easily 
check that the same asymptotic power law f7r(Q^) ^ holds for the %l}{x,k±^ 

functions with any power k, and also for the exponential wave function il}{x., k±) 



x){l + bkl)^] 



e "'"^ I \J x{\ — x). This puzzling 
observation has a rather simple explanation: the valence parton densities in these models with massless quarks are 
constant: /(a;) = 1, and it is this singular j{x)\x^\ const behavior which is responsible for the 1/Q^ contribution 
to the form factor. If we would "correct" the model so that /(x) has a more realistic ~ (1 — x) behavior for x close 
to 1, the Feynman mechanism contribution would have a l/Q'' asymptotic behavior. 

An efficient way to obtain the asymptotic expansion in powers and logarithms of IS? jQ"^ is based on the Mellin 
representation for the denominator factor: 



a: + 7(l-7)(l-xO 



4A2 



1 

2^ 



r(-j)r(j + 3) 



(5— ioo 



{\ - {\ - xY X 



-J-3 



V4A2 



dJ . 



Now, the 7 and x integrals can be calculated in F-functions to give 



j;(Q') = 



''^'^ 6 r(-^)r(i - J) , + 3) V rfJ 

-.-,00 ^J+i)r(2j + 5)^^"^ + '^^^^ + ^^UA^J 



(6.5) 



(6.6) 



The integrand has poles at integer J in the left half-plane. We explicitly displayed the rightmost of these poles l/( J+1). 
It corresponds to x ~ 1 integration and gives the leading asymptotic contribution equal to VlK^ jQ'^ . Expanding the 
integrand in the vicinity of J = —2, —3, eic, we can get subleading contributions. Note that the singularity at J = —2 
is a double pole 1/(J + 2)^. Hence, this contribution will have the {K^ IQ'^Y\vl{Q'^ j t^) term. One can also close the 
integration contour in the right half-plane. This procedure gives the small-Q^ expansion of F^{Q'^), the first terms of 
which are explicitly written in (5.8). 



VII. DOUBLE DISTRIBUTIONS IN TOY MODEL 



Let us now analyze the connection of our expression for the pion form factor 



1 



Ms) Jo ■ Jo 



dx i d-y 



67(1 - 7) 



1 + 



2(£-l/2)f 



>£2(l-x)Qf 



(7.1) 



x(l-x) +7(1 7) 4m^x 

with generalized parton distributions. Introducing the variable y = (1 — a;) 7, we can rewrite this formula as 

6y{l-x-y)/{l-xf 



FAQ^) = 77-^/ dx [ dye{x + y<l 
^\^) Jo Jo 



Xil-X) +2/1-'- X V) ^^■2^(^l_^) 



(7.2) 
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It may be treated as the standard representation [3] 

FAQ^) ^ f dx f F{x,y;-Q^)dy 
Jo Jo 

of the pion form factor in terms of the double distribution^ 

6yil-x-y)/{l~x)^ 



F{x,y;t)^0{x + y< 1) ■ 



(7.3) 



(7.4) 



(we switched to i = —Q^ to conform with the standard notations for generalized parton distributions). This double 
distribution has correct spectral properties [3,5]: it vanishes outside the triangle < x,y,x + y < 1. It also satisfies 
the Munich symmetry condition [16] 



F{x,y;t) = F{x, l-x~y;t) . 
Furthermore, for t = 0, it has the factorized form suggested in Ref. [9] 

F{x, y;t = 0)= e{x + y < 1) h{x, y) f{x) , 
in which the y-dependence appears only in the normalized profile function 

h{x,y)^6yil~x-y)/{l-xf 

satisfying 



h{x,v)dy 1 



The remaining factor 



A{s) 



1 



x{l~x) 



(7.5) 
(7.6) 
(7.7) 

(7.8) 
(7.9) 



depends on x only, and may be interpreted as the parton distribution for the valence quarks inside the pion. For nonzero 
t, the profile function h{x,y) also factorizes out in the expression for the double distribution, Eq.(7.4). However, it is 
multiplied by a function that has a nontrivial dependence on all three variables x, y and t. In Fig. 3, we plot F{x, y; t) 
as a function of x and y for a few values of t. 






FIG. 3. F{x,y;t) as a function of x and y for three values of t = 0, —0.5 and —1 GeV^ 



In Appendix A, it is demonstrated that the variables x, y in Eq. (7.4) have the same meaning as in the DD definition. 
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The form of the double distribution presented above corresponds to parametrization k = xpi + yr of the active 
quark momentum k in terms of the initial pion momentum pi and the momentum transfer r — pi—p2- The final state 
pion has then the momentum P2 = pi — r: the initial and final pions are not treated symmetrically in this formalism. 
To reinforce the symmetry, one should introduce the average momentum P — (pi +p2)/2 (the initial and final pion 
momenta are then P± r/2, see Section II) and write the active quark momentum as = f3P + (1 + a)r/2. To get the 
relevant double distribution a; t), we write the 7 variable as 7 = (1 + ?y)/2 and then introduce a by a = (1 — P)"/]. 
This gives 



f{(3,a;t)^e{\a\<l-(3) 



[(1-/3)' 



']/{l-f3r 



A{s) 



/3(l-/3) 



- [(1 - - a2] 



16m2/3(l-/3) 



(7.10) 



The normalized profile function in this case is |[(1 — /3)^ — a^] /(I — /3)^. 

Integration over k± can be performed in a similar way for a more general power-law LC wave functions ip{x, k±) ~ 
(a + bk'^)~'^ / ^^/x{l — x). In this case, one obtains DDs with K-dependent profiles ~ [(1 — /?)^ — a^]'*~^/(l — 0)'^'^~^ . 
The power of the denominator factor in Eq. (7.10) also changes from 3 to 2k — 1. Note, that the faster the decrease 
of k±) with k±, the narrower is the a-profile of the resulting DD and the faster its decrease with — t. The purely 
exponential wave function (k — s- 00) gives an infinitely narrow profile function 6(a) (or S[y — (x — l)/2] in the case 
of F{x,y;t)). The integral over a (or y) is trivial, and this is the formal reason why the Gaussian model gives a 
one-dimensional integral representation for the form factor. As suspected, the Gaussian model is "too narrow" : it 
cannot reveal the y (a) profile feature inherent to DDs in general case. 



VIII. SKEWED DISTRIBUTIONS IN TOY MODEL 

Having the expression for the double distribution, wc can construct the nonforward distributions [5] in the standard 
way from 

:F^{X,t)^ ( dx ( F{x,y;t)6{X-x-Cy)dy . (8.1) 
Jo Jq 

Note, that for F{x,y;t) given by Eq.(7.4), the integrations again can be performed in elementary functions. Hence, 
in this particular model, one can analytically study the interplay between X, ( and t dependence of the nonforward 
parton distributions (though the expressions are now really lengthy). 




FIG. 4. jr^(X, t) as a function of X and t for three values of ^ = 0.2, 0.4 and 0.6. 



In Fig. 4, we plot Ti^iX, t) as a function of X and t for some values of (■ Note that for each value of C, the nonforward 
distributions satisfy the reduction formula 

/ Tc{X,t)dX ^F^{-t) . (8.2) 
Jo 
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An important point is that the right hand side here has no dependence on the skewedness parameter 
One can also use the symmetric double distribution /(/?, a; t) and the relation 

/■I 

H^^\x,£,;t)= dp r''\(3,a;t)S{x- (3-^a)da (8.3) 

Jo J-l+0 

to obtain Ji's ofF-forward parton distributions Note, that for the infinitely narrow profile function 

f^°°^{l3,a;t) = F{p,t)6{a) corresponding to the purely Gaussian wave function, the OFPD H{x,^;t) is given by 
the ^-independent function H^^'^ {x,^;t) = F{x,t): there are no skewedness effects. Thus, the SPDs obtained from 
the DDs based on power-law Ansatz have a richer structure. 



IX. "REALISTIC" MODEL 



The function f{x), Eq.(7.9) was interpreted above as the toy model version of the valence quark distribution in the 
pion. However, its form strongly differs from the usual phenomenological parametrizations. At a normalization point 
/U ~ 1 GeV, the latter have the form close to 

nx) = l{l-x)/V^, 

with the 1 / a/x reflecting the Regge behavior due to exchanges that are not taken into account in the toy model. In the 
latter, wo assumed that the contribution from the higher Fock components to f{x) has the same shape as the two-body 
one. Also, the expression that wc use for the two-body wave function is again just a model guess. In particular, as 
shown in Appendix A, the double distribution of our toy model can be obtained from the scalar triangle diagram taken 
at spacelike virtualities of the external currents imitating the pions. In the spin-1/2 case, there are extra x-dependent 
factors originating from the numerators of propagators. Thus, one should not take the .x-dependence of the toy model 
DD too seriously. On the other hand, we observed that the ^/-dependence of the model DD has a rather universal 
structure: for t = 0, it is given by the profile function [y{l — x — y)]'^~^ only. For t ^ 0, the y-dependence appears 
also in the denominator factor, but it has a simple structure basically resulting from kinematics. These observations 
suggest to "minimally correct" the model DD: to change its x-shape at t = without changing the pattern of its y 
and t dependence. To preserve the analytic structure of the interplay between the t vs. x and y dependence of DDs 
dictated by the simplest Ansatz (7.4) we take the model 

F^ix,y;t) = F{x,y;t)Q^ . (9.1) 
In terms of the effective two-body-like LC wave function, this corresponds to the change 



V(x, fci) ^ ^""{x, ki_) = i,{x, k^) ^ . (9.2) 

The parameters m and s of the new model should again be fixed by fitting the slope of the pion form factor &t t = 
and its behavior in the — t 1 GeV^ region. In Fig. 5, we show the curve for the pion form factor obtained with the 
Ansatz (9.1) and the values m = 0.46 and s = 0.81. It practically coincides with the curve obtained within the original 
model for m = 0.3 GeV and s = 0.95 in the region of "not-so-high" transfer. 
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0.2 



0.1 



0.2 0.4 0.6 0.8 1 2 4 6 8 10 

FIG. 5. Form factor -F,r(Q^) (left) and Q^F^((5^) GeV^ (right) obtained witii the "realistic" model for double distributions 
(solid line). For comparison, we present the results for the original model (dashed) and "p - meson fit" (5^/(1 + Q^/(0. 77 GeV)^) 
(dotted). 

With the new DDs, one can obtain a realistic model for SPDs J-^{X,t) via (8.1). The SPD is presented in Fig. 6 
as a function of X and t for some values of C- 




FIG. 6. J-^{X,t) obtained with the "realistic" DD as a function of X and t for three values of = 0.2, 0.4 and 0.6. 



As a more explicit illustration of the i-dependence of SPDs, in Fig. 7 we show SPDs F(^{X; t) with different C's for 
two different values of both for the original and the "realistic" model. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



FIG. 7. SPDs Fc^{X;t) with C, = 0.1, 0.2, 0.4, 0.6 are shown for two different values t = (solid lines) and t = -0.2 GeV^, 
both for the original (left) and "realistic" (right) model. The curves corresponding to larger C, have maxima at higher X. 

When C, increases, the maxima of SPDs shift to higher values of X . The rate of change is more drastic in the case of 
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the realistic model, where SPD changes from a monotonically decreasing curve for C = (corresponding to the usual 
parton density) to a shape resembling that of distribution amplitudes, as ^ tends to 1. It is interesting to analyze the 
limiting case C 1- As demonstrated by M. Polyakov [23], in the soft pion limit, — > 0, C ^ 17^ = 0, the isovector 
part of the pion SPD should coincide with the pion distribution amplitude. To check if this constraint is satisfied by 
our models, we take the function Fi^=i{X;t — 0) and symmetrize it with respect to X 1 — X to project onto the 
isovector component. In Fig. 8, we show the results both for the toy model and the realistic model. For the toy model, 
we obtain a double humped curve, resembhng the Chernyak-Zhitnitsky (CZ) model (p^^{X) = 30 (1 - 2X)^X{1 - X). 
More precisely, this curve can be fit, with a good accuracy, by the sum 0A3 Lp°'''{X) + 0.57 ip'~"'^ {X) of the CZ and the 
asymptotic distribution amplitude Lp°-'^{X) = 6X{1 — X). The pion distribution amplitude in the toy model can be 
obtained directly by integrating the two-body wave function ip(x,k±) over the transverse momentum. The result is close 
to the asymptotic amplitude, so one can say that the toy model does not satisfy the constraint imposed by the Polyakov 
soft pion theorem. On the other hand, in the realistic model, the function [FQ=i{X;t = 0) + ~ X;t — 0)]/2 

is very close to the asymptotic form and to the distribution amplitude obtained from the two-body wave function. 
Experimentally, the pion DA is known to be rather close to the asymptotic form. Thus, the realistic model de facto 
satisfies the constraint imposed by the soft pion theorem. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



FIG. 8. Left: isovector part of F(^i{X;t = 0) for the toy model (solid) compared to the asymptotic pion DA (dash-dotted), 
CZ distribution amplitude (dash-double-dotted) and DA obtained from the two-body wave function (dashed). Right: isovector 
part of F(^i{X;t = 0) for the realistic model (solid) compared to the asymptotic pion DA (dash-dotted) and DA obtained from 
two-body wave function (dashed). 

X. CONCLUSIONS 

In this paper, we demonstrated how to obtain a model for the valence (or C-odd) pion double distribution F{x, y; t) 
and the skewed parton distribution J-Q{X;t) satisfying, by construction, such important constraints as reduction 
relations to usual parton densities and form factors, spectral and polynomiality conditions. SPDs derived in our 
model have a nontrivial interplay between X, C, and t dependence. Furthermore, they were adjusted to describe pion 
form factor for all available t, so we expect that our model describes the t-dependence of the pion GPDs of the valence 
quarks both for small and large t. The ability to have a unified model for GPDs from t = to |i| ~ lOGeV^ is 
especially important in (future) applications to nucleons, for which GPDs are already being studied experimentally 
both for small (DVCS) and large t (wide-angle Compton scattering). 
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APPENDIX A: DD IN A SCALAR MODEL 



Our conversion of the integral representation for the form factor F(t) into a function F{x, y; t) of three variables 
may look like a rather ambiguous exercize. Below, by a covariant field-theoretic calculation, we demonstrate that x 
and y really have the meaning of the variables of a double distribution. 



FIG. 9. Box diagram for DVCS in a scalar model. 



First, consider a one-loop box diagram for the scalar analog of deeply virtual Compton scattering amplitude (see 
Fig. 9). The initial and final "pions" are imitated by scalar currents 11 corresponding to spacelike momenta pi and 
P2, the initial "photon" momentum is q\, and that of the final one is q2. The momentum invariants describing this 
4-point function are 



2 2 
-9l,92 



0, s = (pi + qif, t={pi- p2f 



Using the a-representation 



le 



'f 

Jo 



'+^^) da 



(Al) 



(A2) 



for each of four scalar propagators and calculating the resulting Gaussian integral over the loop momentum k we 
obtain 



T{pi,P2,qi) = 



exp 



Q!i + a2 + as + a4 



p (m^ — ie) 



da 1 da2 da^da^ 



(A3) 



where p = ai + a.2 + a:^ + a^. Wc arc interested in the Bjorken kinematics when there are two large variables s 
and which have the same order of magnitude s ~ (l/xsj — while other invariants are small. The large-Q^ 

asymptotics in this situation is determined by integration over the region where the coefficients accompanying s and 
vanish simultaneously. Otherwise, the integrand rapidly oscillates and the result of integration is exponentially 
suppressed. Integration over ai ~ region is the simplest (and, as inspection shows, the leading) possibility. It 
corresponds to hard momentum flow through the propagator connecting the photon vertices. Performing the ai 
integration, we obtain 



T{pi,P2,qi) 



(•i(i2(/fi;jr/n4/' A" 
saz/X — Q^a^/X + ie 



exp <^ - [a2a4t + aaiaApj + a2pl)] - iX{m^ + 0{1/Q'^) , (A4) 
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where A = a2 + 0:3 + 04. Denoting v = 2{piqi) and writing s = ly—Q^ {pf is neglected compared to Q^), we represent 
the Q^-dependent term in the denominator as va^/X — Q^{as + a4)/A, or finally as ulas/X — ({1 — a2/X)], where 
C = jv- Introducing the double distribution 



P{x,y;t,pl,pl) = i J 6 (x- (5 - ^) |x ["20:4^ + Q;3(q;4P? + "2^2)] - iX{m^ - ie) 



A2 

(A5) 

we can write the scalar DVCS amplitude in the partonic form 



= -- f f 

^ Jo Jo 



rri \ -II P{x,y;t,pl,pl) 

T{pi,P2,qi) = — / / — T — j—^dxdy. (A6) 

x + yC-C + ie 



A few comments are in order. First, we reemphasize that there is no C-dependence in the definition of DDs. The 

second comment concerns the spectral properties of DDs. It is easy to see that both variables x, y vary between and 
1. Furthermore, their sum is also confined within these limits: < x + y < 1. Finally, the hard amplitude depends on 
the DD variables x, y through the combination x + y( only, so one can treat it as a new variable X = x + y( and use 
nonforward parton distributions ^q{X; t) instead of DDs F{x, y; t). 

With the same technique, one can calculate the nonforward matrix element {p2\On\pi) of the composite operator 
On = ip{izdy^f and obtain its a-representation 



Jo 



fa3{piz) + a2{rz)Y hr ., , 2, 2m -^i 1 . .\ da2da3da4, 

I I e^p<j[a2ait + a3{aipt + a2pi)\-iX{m^-ie)> • l^^) 



Note, that that the derivative {izd) acting on field (p is expected to give the factor (kz), where k is the momentum 
of ip. Eq.(A7) shows that (kz) = a3{piz) / X + a2{rz)/X, i.e., a^/X and a2/A should be interpreted as the variables x 
and y of the double distribution F{x,y;t). Alternatively, using the binomial expansion for the (. . factor, one can 
sec that the coefficients in front of (pi2;)"~'^(rz)'' in this expansion are given by the x'^~^y^ moments of the DD (A5). 
Integrating over A, we get the relevant DD explicitly 

F{x,y;t,pl,pl) = ~{y{l-x-y)t + x[{l-x-y)pl+ypl]-m']~^ . (A8) 

Putting the "pions" on equal footing by setting p\= p'2 = —M^, we get a DD 

F'^^\x,y;t,M'^) = {-y{l - x - y)t + x{l - x)M'^ + m'^}~'^ (A9) 

satisfying the y ^ 1 — x — y Munich symmetry condition. This expression can be rewritten in the form 

resembling the DDs obtained from the power-law wave function. Introducing A-^ = M^/4-|-m^, we get the expression 

F(^)(x,,;,M2)U..4(A^_.) = ^ |l + J-^^J^^^ - (All) 

whose denominator factor has the structure close to that of the DDs obtained in the model with the power-law wave 
functions. However, the denominator power is (—1) instead of (—3). Applying {p\d / dp\){p2d / dp^) to F{x, y; i,Pi,pi), 
Eq. (AS)*, and setting pl=pl = -M^, we obtain 



tit is easy to see that the IlipLp scalar vertex differentiated with respect to the virtuality of the scalar H-current corresponds 
to the K = 2 power-law wave function (4.4). 
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Now, using = + m^, we end up with the DD differing from the toy model DD, Eq. (7.4) just by the 

a;-dependent factor 1/x and an overall normalization. Note, that our toy DD was based on the formula for the vector 

form factor, hence, for full correspondence, we should consider DD related to operators containing the extra i 
derivative. This results in the extra factor of 

y (K + ^^y^ K - P^) = 22;^" + {2y-x- l)r'* . 

As expected, the P'^-part contains the missing factor of x. Since the r'^-part is Munich-antisymmetric, it does not 
contribute to form factor and forward densities (sec also [24]). In general, such terms arc not restricted by the 
reduction relations (2.13), (2.14). However, they contribute to SPDs for C 0, and their modeling deserves a separate 
consideration. 



APPENDIX B: SPDS IN THE IMPACT PARAMETER REPRESENTATION 

In this Appendix, we investigate the properties of our model SPDs in the impact parameter representation. For the 
scalar triangle diagram, such an analysis was recently performed by Pobylitsa [25]. He also uses the a-representation 
and double distributions for the triangle diagram, but takes the version of double distributions in which the plus 
component of the momentum of the spectator system is written as upf + ■ Instead of Eq.(A5), we would have 
then 



P{u,v;t,pl,pl) = i J (" ^ y) ^ y) cxpl^j[a2a4t + a3{a4p'i + a2pl)] -iXim^ -ie)^ 



da2dazdai 



A2 

(Bl) 

This DD is related to SPD H{x,i:t) by 

H{x,^;t,pl,pl)= f f S{l-x-u{l + ^)-v{l-^))P{u,v;t,pl,pl)e{0<u + v<l)dudv . (B2) 
Jo Jo 

Now, one should take Pi = P2 = to^> * = — (|A_l|^ + 4^^m^)/(l — ^^) and calculate the double Fourier transform 

B(«;,J./||i,(^....;,(,,,_^4±|M) ^ (B3, 

The i5-function in Eq.(B2) can be rewritten as S{1 — u/ri — ?;/r2)/(l — i), where the parameters ri, r2 given by 
ri = (1 — 5:)/(l + ^), r2 = {1 — x)/{l — ^) have the meaning of the spectator's plus-momentum measured in units of 
the initial or final pion plus-momenta. Due to this ^-function, we can write u = zr\,v = (1 — z)r2, with < 2; < 1 in 
the X > ^ region. Finally, the integral over dXdz can be transformed into integration over the variables di = zX and 
0-2 = (1 — z)X. A remarkable fact is that the resulting integrand /(cri,cr2) factorizes I{a\,a2) = Ji{cri)J2{cr2)- As a 
consequence, the expression for B{x, ^; b_\_) also has a factorized form [25] 

B{x,^;b±) = ^Vo{n,{l-Ob±)Vo{r2,{l+Ob±) , (B4) 



where the generalized impact-parameter LC wave functions 

a{m^ — r(l — r)mj) 



I r^da 

Vb(r,cx) = -; — / — exp 
A-KT Jo a 



znr \ r 



m2 - r(l - r)ml) (B5) 
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can be expressed through the modified Bessel function Kq. As demonstrated in Refs. [25], the factorized representation 
(B4) guarantees that the positivity bounds [26,9,27,14,28-31] for this SPD are satisfied in the model with scalar quarks. 
The same SPD multiplied by (1 — x) satisfies the positivity bounds for spin-1/2 quarks. 

If we would proceed as in Appendix A, i.e., first differentiate Eq.(Bl) with respect to Pi and P2 and then take 
Pi = pI — ^^Jj w6 would get an extra factor Q!2Q!4a§/A^. After the transformations described above, this would result 
in the factor rir2C7iC72[l — (ricri + r2<J2)/{'^\ + 0'2)]^, and the integral over (Ji,(J2 cannot be factorized into a product 
of two separate integrals over cri , (T2 • The unfactorizablc piece comes from the a\ factor resulting from differentiation 
with respect to external virtualities. To avoid this factor, but still preserve the ~ 1/^1 behavior of the effective 
IMF wave function, one can perform differentiation with respect to the squares of the quark masses {i.e., take all 
the masses different, differentiate with respect to 's corresponding to lines 2 and 4 and then take all the masses 
equal). This would produce the factor a2ai, or eventually rir20iG2, which does not violate the factorized structure 
of the integrand. In the impact parameter representation, the result has the structure of Eq.(B4), but with Vo(r, c±_) 
substituted by the expression 

V,{r,c^) = ====Ki (\^^i^r{l-r)ml/mA (B6) 

involving the modified Bessel function K\. For the original IMF wave function, differentiation with respect to the 
active quark mass is equivalent to choosing the Ansatz 

k±) = , ,^ (B7) 

.'Ei/.'j;(l — xj[a + bkj_\'^ 

instead of Eq.(4.4). It has the extra 1/x factor enhancing the wave function at small x. In the spirit of our discussion 
of the a;-dependence in the main text of the paper, we may say that such a function more adequately models the 
contribution of higher Fock components. 

The positivity bounds are satisfied also in a more general case when B{x, ^; b±) is given by a sum [32,25] 

B{x, ^; b^) = (1 - xf+' Qnin, (1 + ^)b±)Q*n{r2, (1 - Ob±) , (B8) 

n 

where N = ior "scalar quarks" and = 1 for spin-1/2 case. This opens a possibility to build models for GPDs con- 
sistent both with the polynomiality and positivity constraints. The simplest idea is to start with the a-representation 
(A5) for DD corresponding to the scalar triangle diagram, and modify it by multiplying the integrand by a fmiction 
R{w?y/a2a4) depending only on the product a20n (choosing the argument as -^0204 we get eventually a function 
of {1 — x); the parameter m? was included to make the argument of R{a) dimensionless). If this function can be 
represented as 

R[a) = ^ i?„o" (B9) 

n 

with all Rn positive (the sum should be imdcrstood in a wide sense, it can involve integrations), the function B{x, ^; &j^) 
in such a model would have the structure of Eq.(B8), and positivity constraints are satisfied. The model double 
distribution based on R{a) gives the following expression for H{x,£,;t) in the x > region 

H{x,^;t)\x>^ = {1 - x)'^~^rir2 dp / dz R{p^/zzrir2) 

Jo Jo 

X exp {— p [1 — (1 — zri — zr2){zri + zr2)m'^/m^ — zzrir2t/m^]} (BIO) 

(we use the notation z = 1 — z here, and later we also use x = 1 — x, etc.). Note that we performed Wick rotation 
aj —iaj in the original Q-rcpresentation (Bl), which is justified if the pion stability condition < 4m^ is satisfied. 
In the forward limit = 0,t = 0) this gives 

pi POO 

f{x) = {l-x)^+^ dz rfpi?(p\/^S)e-''(i-^^'"'/'"') . (BU) 
Jo Jo 



20 



The function R{a) should be adjusted to fit experimental forward distribution /^^^{x), and then it can be used for 
calculation of H{x,^;t). In particular, for massless pion, the coefficients Rn can be expressed directly 

n + 1 

[r(n/2 + l)P^"+^+^ 
in terms of the coefficients of the x'' expansion of the forward distribution 

r^{x) = Y,M^-x)'' . (B13) 

k 

Note, that for the simple model f^{x) = |(l-a;)/v^ of Section IX, all the coefficients Ak = r(A:- l/2)/r(l/2)(A;- 1)! 
are positive. Alternatively, since both Eqs.(BlO) and (Bll) in the t = 0, = limit involve the same integral of the 
i?- function, the only change being x — > y^rlri, in this case we can directly write H through f{x) 

iJ(i,?;i = o)U>i;™j=o=(^)'^ '/(i-VfI?^) = (i-a^^-'^/v(|i-^^==) • (B14) 

In the case of nonforward distributions, we have 

J-(X,C;t = 0)U>C;r»|=o = (l-C)^^-^)/'/(l-^^) • (B15) 

Taking f{x) = f^{x) and N = 1 (spinor quarks), we get the curves which are very close to the X > parts of 
the realistic model curves shown in Fig. 7. In the region X < (, the functions can be obtained only from formulas 
explicitly involving double distributions. In particular, the scalar R{a) model gives 



POO 

F{x, y;t)= / e-^f R[p ^y{l-x-y)] dp , 
Jo 



(B16) 



where p{x,y,t) = 1 — y{l — x — y)t/m? — x{l — x)m'^/m^. Fixing R{a), e.g., by the requirement that /(x) = f^{x) 
in the — > limit (which allows to get the result in analytic form) and using Eq.(B12) we get 

Fix, y,t)-^ A„+, r2(„/2 + i)[p(,^^^,)]„+i • (B17) 
This expression can be also written as 

-1 ra+l 



F{x,y;t) = Y,An+ih^'''^\x,y) 



(1 



.P{x,y,t) 



(B18) 



where h^'^^'^\x,y) is the normalized profile function of order n/2. Thus, DD F{x,y:t) in this model is given by a 
sum of powerlike terms similar to those discussed in Section VII. The main difference is that the y-profile now is 
not universal: one should expand the forward distribution f{x) into a power series over (1 — ./ ) and supplement the 
(1 — a;)""*"^ term by the n-dependent profile function /i^"/^^ (x, y). Because of the correlation between the power of 
(1 — x) and the order of the profile function, the shape of the y-profile of the double distribution changes with x. Since 
all parton distributions f{x) tend to infinity as x goes to 0, the small- a; region is dominated by terms with large n, and 
the profile is more narrow when a; ^ 0, becoming infinitely narrow as a; ^ 0. Note, that in the case when the proffie 
is infinitely narrow for all x, i.e., when F{x,y,t = 0) = S[y — (1 — x)/2\f{x), the t = skewed distribution is given 
by ^^{X) = f[{X — (/2)/{l — (^/2)]/vT^-C'- it repeats the form of the forward distibution f{x) and is infinite for 
X = (/2. To check if this also happens for the model (B18) with the changing profile and realistic ~ x~°'^ behavior for 
small X, we took /°'^p(x) = f^{x) (or, what is the same, A„+i = r(n + l/2)/r(l/2) n! in Eq.(B17)) and constructed 
model skewed distributions ^({X) both in X > ^ and X < ( regions (see Fig. 10). The resulting functions are indeed 
singular for X = (/2, but finite otherwise. 
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FIG. 10. SPDs t = 0) with C = 0.2, 0.4, 0.6 obtained from Eqs.(B18) and (B13). The forward distribution was modeled 

by /■^(a;) = 1(1 - x)/y^. The curves tend to oo for X = C/2. 

We plan to perform a more detailed study of GPDs within the i?(a)-model in a separate paper. 

APPENDIX C: TWO-PION DISTRIBUTION AMPLITUDE IN SCALAR MODEL 

Another referee's request was to apply our approach to processes of two-pion production in 7*7 — > tttt or 77 tttt 
collisions. The 7*7 tttt process in the kinematics when 7* is highly virtual, while s = rri^^ is small [33] is the crossed 
channel reaction to DVCS while 77 tttt process in the kinematics when s, \t\, \u\ are large [34] is the crossed channel 
reaction to wide-angle Compton scattering. Here we are going to consider only the simpler case of the kinematics of 
the first process. 



- VP, 



(l-z)P \P,= (1-0P 

FIG. 11. Plus-momentum flux structure of the two-pion distribution amplitude. 

Originally [33], it was proposed to describe the nonperturbative stage of this process by two-pion distribution 
amplitude (27rDA) $(2;, C; s) which describes the conversion of two quarks with plus-momenta zP+ and zP^ into 
two pions with momenta = C^*^ ^nd — ^P~^ (see Fig. 11), where P is the total momentum of the pion pair 
(recall that the invariant mass of the pair is small: = s <C IGeV^). Later, Teryaev [35] proposed to use the 
double distribution description (see also Ref. [36] for further developments) corresponding to parametrization of the 
plus-component of the spectator momentum as upf — (note, that both momenta are now outgoing, and this is 
reflected in the change of the relative sign of pi and p2 parts of the spectator momentum compared to the DVCS 
case). In P+ units, the spectator plus-momentum can be written either as (C — z)P+ (using the 27rDA variable z) or 
as u(^P~^ — v^P^ (using the DD variables u, v). Thus, the connection between the two sets of variables is given by 

z — {1 ~ u ~ f )C -I- V . 
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So, we can express the 27rDA ^{z, s) in terms of the DD M{u, v; s): 

^z,C;s)= f du f dve{u + V < l)d{z - C{1 - u- v) - v) M{u,v;s) . (CI) 

Jo JQ 

The DD representation for ^{z, (; s) can be used to derive some general properties of 27rDAs, such as the polynomiahty 
condition 



/' 

^0 



^"$(^,C;s)rf^ = ^ i^iC' , (C2) 



1=0 



which states that z" moment of ^{z, (; s) is nth order polynomial of (. 

To model 27rDAs by superposition of perturbative contributions, we start with the a-representation for the relevant 
scalar diagram. The double distribution can be written similarly to Eq.(Bl) 

M (w, v;s)=im J o — —j o — —j exp | - [a2a4,s + as^ai + a2jm^\ - iX{m - le) > . (C3) 

Here, from the beginning we take Pi=P2 = m^, and add the overall factor to make the function dimensionless. 
Note, that for the triangle diagram we have {1 — u — v) = as/ X = ps. So, we write the scalar triangle version of the 
27rDA as 

poo pl pl 

^z, C; s) = im^ / dX / dp2 / d/Ss e{p2 + p3<l) S{z - C/?3 - ^2) 
Jo Jo Jo 

exp {iX [(1 -02- /?3)/?2S + /33(1 - 03)ml -m^ + ie]} . (C4) 
Integrating over A and incorporating the (5-function to calculate the P2 integral, we get 

/.min{0/C,z/C} _ 

Hz, C; s) = / dps [l-{z- PsOiz - PsOs/m^ - Psil - PsWJm^ - ie] . (C5) 

This representation explicitly demonstrates the well-known fact (see, e.g., [37,36]) that $(z, C,; s) is non-analytic at the 
point z = Q. The integral can be taken in general case, but it is instructive to analyze the simplest limit s = , m\ = 0. 
In this case, the result is the function 

Hz, C; s = Q)mi=o = J ^(^ < + 0{z > (C6) 

that coincides with a part of the pion DA evolution kernel. Its eigenfunctions are the Gegenbauer polynomials 
Cn'^{2( — 1) and the eigenvalues are l/(n -|- l)(n -|- 2) [38,39]. Hence, we can write 

Hz, C; s = 0)^. .0 = 4z(l ~z)J2 fa + 1)2! + 2)^ " ^^^"^'(^C - 1) • (C7) 

n=0 ^ ^ ^ ' 

It is convenient to write 27rDA as a sum over z{\-z)cT{^z-\), since these are the eigenfunctions of the evolution 
kernel. On the other hand, the combination (2C — 1) is related to the cosine of the angle between the pions' momenta, 
so it is natural to expand the ^-dependence of <i>(z, C; s) in the Legendre polynomials Pii^C, — 1). Using the formula 

c\i\x)-c\^^{x) = {:ii^m{x), 

we can write C^/'(2C-1) as a sum of P;(2C — 1) and obtain 

Hz, C; s = o)„. =0 = 4.(1 -z)Y: (^ + iJ^ + 2)^ ^"^'(2" - ^) S + - ^) 2 • 

n=0 V V ^ (=0 
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This expansion has the structure of the general representation for 27rDAs proposed by Polyakov [23]. In specific 
models, only the first terms of the expansion are included. It is easy to check that, in our case, the exact result (C6) is 
well approximated by a few first terms of the Gegenbauer expansion (C7), even in the vicinity of the non-analyticity 
point z = (. Polyakov [23] considers also a more complicated s ^ case, using the tttt scattering information to model 
the s-dependence. Recently, Kivel and Polyakov [36] used chiral perturbation theory to include 0(mj) corrections to 
the chiral limit. 

Now we want to show how one can incorporate information about usual (forward) parton densities to build models 
for 27rDAs. To this end, it is convenient to write the pion momenta as pi = P/2 + r and p2 = P/2 — r (the plus- 
components are implied, but we omit the + superscript here and below). The quark momenta can be written then 
as 

ki = ^-^P + xr , k2= ^—^P - xr , (C9) 

where the variables a and x are related to u,v hy x = 1 — u — v and a = v — u. The support region is jaj < 1 — x. 
The 27rDA $(2:, C,; s) is related to the double distribution F{x, a; s) by 



^z,C;s)= [ dx [ daS{z - 1/2- x{C- 1/2)- a/2) F{x,a;s) . (CIO) 

Jo J-l+x 

In this description, the total pair momentum P is shared by the quarks in fractions (1 + a)/2 and (1 — a)/ 2, while 
the relative momentum r is carried by active quarks in fractions x and —x. Hence, the relevant double distribution 
F(x, a; s) is the timclikc analog of the function f{x, a; t) considered in Section II. In the forward limit, it reduces to 
the usual parton densities: 



l-x 

F{x,a;s = 0)da = f{x) . (Cll) 

l+x 



In the = case, we have F{x, a;s = 0) = -^O^al < 1 — x), hence the integral in (Cll) gives (1 — .t), which is exactly 
the forward distribution for the scalar massless triangle. To get a more realistic f{x), we can use the i?(a)-model 
described in Appendix B. It gives 

pinin{z/(,z/Q POO f r _ 1 1 / I ^\ 

$«(^, C; s) = y dx J exp |-p ^1 - (^ - xC){z -xC)^-xx^^^ R [p {z - xC){z - xQ j dp . (C12) 

According to Appendix B, the i?(a)-model is equivalent to a sum of "wave function overlap" contributions of Eq.(B8) 
type, similar to those obtained within the light-cone approaches (see, e.g., [40,31,41]). However, the R{a) construction 
has the advantage that it provides also a model for 27rDA. In the standard light-cone formalisms, the 27rDA would 
involve the q ^ nq vertices which cannot be interpreted as light-cone wave functions. 

In the scalar triangle model, the 27rDA ^{z, C; s = 0) is obtained from the same DD F{x, a;s = 0) = f{x, a;t = 0) 
which produces the t = OFPD H{x, ^; t = 0). Comparing Eqs.(ClO) and (8.3), we can formally write 

<p{z,C,s = 0)=H(^^,-J--;t = o]. (C13) 



.2C- 1 2C 
The relation is even simpler 

^{x, i;s = 0) = H (x/i, 1/i; t = O) (C14) 

for the 27rDA (j){x,^;s = 0) written in the symmetric variables x = 2z — 1 and ^ = 2( — 1. Since |^| < 1, the OFPD 
H{x/^ , 1/^; t = 0) is taken at the skewedness values with absolute magnitude larger than 1. Hence, as suggested by 
Teryaev [35], the 27rDA may be treated as a continuation of OFPD into the |^| > 1 region. More precisely, -ff|t=o and 

$|s=o may be treated as j^j < 1 and j^j > 1 components of the same function. 

To make parallel with / = and 1=1 components of 27rDAs in QCD, one should take combinations 
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(C15) 



which are symmetric or antisymmetric with respect to the middle point z = 1/2. They are given by summation over 
even or odd n in Eq.(C8). Taking s — 0, = and fixing R{a) in the same way as in Appendix B, we obtained the 
curves shown in Fig. 12 . 
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FIG. 12. Two-pion distribution amplitudes ^'^{z,(;s = 0) and $ (z,i^;s = 0) with ( — 0.1, 0.2, 0.4 obtained in the scalar 
7?(a)-model. The forward distribution was modeled by f^{x) = |(1 - x)/y^. 

At 2 = 1/2, the symmetric function in this model is infinite. This result is similar to singularity of SPDs J-(^{X) 
for X — C/2 observed in Appendix B. It reflects the fact that the profile of F{x, a,s — 0) in the i?(a)-model becomes 
infinitely narrow as a; — *■ 0. Indeed, for DDs with infinitely narrow profile for all x, i.e., for F{x, a, s = 0) = f{x)S{a), 
we would have s = 0) = f[{z - l/2)/(C - l/2)]/(l - C/2), which gives infinite result for z = 1/2 if /(O) oo . 

The curves also have cusps for z = C and z — 1 ~ They appear because the DD F{x, a,s = 0) of the i?(a)-model 
does not vanish at the upper corner a; = 0, a = 1 of the support region. This is because the profile function for the 
lowest term of the i?(a) expansion is h'^'^\x,y) = 1/(1 — x): unlike profile functions h''^^{x,y) with n > 0, it does not 
vanish at the border lines x + \a\ — 1. 

Note that 27rDAs of the purely scalar model are symmetric with respect to the change {z ^ 1 — z,C, 1~C} while 
in QCD the 27rDAs describing transition of spin-1/2 quarks into pions changes sign after this transformation. The 
triangle perturbative contributions for this case were considered a few years ago by Polyakov and Weiss [37]. We plan 
to extend their calculation by combining it with the ideas of the present paper. 
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